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Abstract We review recently developed models of galactic discrete sources 
of high energy neutrinos. Some of them are based on a simple rescaling of 
the TeV 7-ray fluxes from recently detected galactic sources, such as, shell- 
type supernova remnants or pulsar wind nebulae. Others present detailed 
and originally performed modeling of processes occurring close to compact 
objects, i.e. neutron stars and low mass black holes, which are supposed 
to accelerate hadrons close to dense matter and radiation fields. Most of 
the models considered in this review optimistically assume that the energy 
content in relativistic hadrons is equal to a significant part of the maximum 
observable power output in specific sources, i.e. typically ~ 10%. This may 
give a large overestimation of the neutrino fluxes. This is the case of models 
which postulate neutrino production in hadron-photon collisions already at 
the acceleration place, due to the likely pair plasma domination. Mod- 
els postulating neutrino production in hadron-hadron collisions avoid such 
problems and therefore seem to be more promising. The neutrino telescopes 
currently taking data have not detected any excess from discrete sources 
yet, although some models could already be constrained by the limits they 
are providing. 
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1 Introduction 

Galactic sources can potentially produce interesting event rates in neutrino 
telescopes. Since they are at shorter distances to the Earth (~ 1 10 kpc) 
compared to extra-galactic sources, the source luminosity required for a 
galactic source to produce the same event rate as an extra-galactic one, is 
orders of magnitude smaller (see Sec. [2|l - Given the large photon luminosities 
observed from some of the sources, it is possible to single out interesting 
candidate neutrino emitters in the Galaxy. A rough estimate of the source 
luminosity required to produce a certain event rate in a neutrino telescope 
(Halzcn 2003) is shown in Sec. 

In order to produce high energy neutrino fluxes, galactic sources must 
accelerate particles at sufficiently high energies. Hillas derived the maximum 
energy E at which a particle of charge Z can be accelerated, from the simple 
argument that the Larmor radius of the particle should be smaller than the 
size of the acceleration region (Hillas 1984). If energy losses inside sources 
are neglected, this maximum energy E (in units of 10 18 eV) is related to 
the strength of the magnetic field B (in units of /iGauss) and the size of the 
accelerating region R (in units of kpc) by the following relationship: 

Ei 8 ~ j3 Z B^q R-kpc (1) 

where f3 is the velocity of the shock wave or the acceleration mechanism 
efficiency. Hence, the maximum energy up to which particles can be accel- 
erated depends on the BR product. Particle acceleration may occur in many 
candidate sites, with sizes ranging from kilometers to megaparsecs. 

Most of the models discussed below assume that neutrinos are produced 
by charged pion decays which in turn are produced in hadron-hadron (pp) 
and hadron-photon (p7) collisions. As a matter of fact, in order to produce 
neutrinos with energies of the order of 10 TeV, the primary hadron energy 
has to be about two orders of magnitude larger in the case of pp collisions, 
and about one order of magnitude larger in the case of P7 collisions. Neu- 
trinos produced as decay products of neutrons (which originate in hadronic 
collisions) are usually much less important since they are produced with 
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much lower multiplicities (with respect to charged pions from pp collisions) 
and typically take only ~ 10~ 3 of the initial neutron energy. Therefore, the 
fluxes of neutrinos from neutron decay are significantly lower than those pro- 
duced in hadronic collisions (through meson decays) due to the steepness 
of the spectrum of accelerated hadrons (described usually by a E~ 2 power 
law). Only in some special conditions interesting fluxes of neutrinos from 
decay of neutrons might be expected, e.g. in the case of neutrons produced 
in heavy nuclei photo-desintegration which do not lose efficiently energy in 
pion production in the radiation field. 

In this review paper, we concentrate on the models of galactic neutrino 
sources which have been developed over the past few years. For reviews on 
earlier works, the reader is referred to, e.g., Gaisser et al. (1995) and Learned 
& Mannheim (2000). In Section |3 we consider high energy neutrino fluxes 
produced during the early phase of supernova explosions, whereas Section^] 
will be devoted to pulsar wind nebulae. Shell-type supernova remnants are 
also good candidate sources for high energy neutrino emission, and this will 
be illustrated in Section [3] In Section [B] we will discuss neutrino emission 
from pulsars close to high density regions, such as those observed in the 
Galactic Center or massive stellar associations, and from neutron stars in 
binary systems in Section Finally, in Sections |H1 and El we will discuss 
models of microquasars and magnetars. Section ^| will be devoted to our 
conclusions. 

The predicted event rates by the models reviewed in this paper should 
vary in the presence of neutrino oscillations. The expected flavor ratios at 
sources of high energy cosmic neutrino fluxes from pp or pry collisions are 
4> Va : <j) v : 4>v T = 1 : 2 :< 10~ 5 . Given the solar neutrino results (see 
for instance Fogli et al. 2004), and particularly the recent results by SNO 
(Ahmed et al. 2004) with best fit oscillation parameters Am 2 = 7.1 • 10~ 5 
eV 2 and 8 — 32.5, the strong evidence for and v r maximal mixing from 
atmospheric neutrino experiments which prefer Am 2 ~ 2 • 10~ 3 eV 2 (Supcr- 
Kamiokande (Fukuda et al. 1998), MACRO (Ambrosio et al. 1998) and 
Soudan 2 (Sanchez et al. 2003)), and the limit on the mixing matrix element 
\U e s\ <C 1 from CHOOZ (Apollonio et al. 2003), these flavor ratios are 
expected to transform into 1:1:1 (Learned & Pakvasa 1995, Athar et 
al. 2000) along > kpc distances in a direct hierarchy scenario. For models 
considering electron antineutrinos from neutron decay (such as Anchordoqui 
et al. 2003a) described in Sec. lfi.2|l . about 60% of the electron antineutrinos 
should reach the Earth and ~ 40% should oscillate into and v T in about 
the same percentage. 

In this paper we give an estimate of the muon neutrino event rates in 
the absence of neutrino oscillations. If neutrino oscillations are to be con- 
sidered, a part of the muon neutrinos from mesons decay convert into other 
flavors and could be detected. Nevertheless, the estimate of event rates for 
t events in neutrino telescopes needs a calculation including the regenera- 
tion effect in the Earth due to the fast tau decay (see Bugaev et al. 2004 
and references therein). Moreover, it should be considered that the topolo- 
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gies induced by v r and v e arc different than induced events and hence 
have different detection efficiencies. As a matter of fact, neutrino telescopes 
measure neutrinos indirectly thanks to their interactions with the matter 
inside or surrounding instrumented regions (for a review on recent results 
sec Montaruli 2003). The new generation of neutrino telescopes aims at the 
detection of the Cherenkov light produced by charged secondaries induced 
in neutrino interactions. 3-D arrays of optical modules, pressure resistant 
glass spheres containing phototubes (PMTs), are deployed in ice or sea/lake 
water depths. For muon neutrinos (and also for extremely high energy v T ), 
the effective target mass is much larger than the instrumented region, since 
the technique profits of the increase of the muon (tau) range with energy. On 
the other hand, neutrinos of other flavors and neutral current interactions 
produce showers in the instrumented regions. 



2 Event rates and luminosities from galactic sources 

For an isotropic emission from a source at a distance D, the relation between 
the luminosity and the neutrino energy flux /„ is: 

L v = AttD 2 J v . (2) 

A source of energy flux /„ in neutrinos of energy E v will produce N neu- 
trino induced upward-going muons in a detector with effective area A and 
exposure time T: 

jv = Ap ^r, (3) 

where the probability P v ^p. that a neutrino produces a muon above a certain 
threshold depends on the neutrino charged current (CC) interaction cross 
section, the muon range and should account for the Earth shadowing effect 
1 . At 100 TeV the probability that a neutrino produces a muon of energy 
larger than 10 GeV is P v ^^ ~ 10~ 4 . Hence, the neutrino energy flux and 
the luminosity required to produce N detectable events are 

and 

^ 6 . 10 ^(£)\(^)(J_)-'e rgs -.. (5, 

This means that in one year the detection of 10 events of 100 TeV neutrinos 
coming from a source at D = 5 kpc in a 0.1 km 2 volume, requires a lumi- 
nosity of the order of L v <~ 10 36 erg/s. This large value of the luminosity 
strongly limits galactic source candidates. 

1 Earth shadowing effects are small for galactic models that predict neutrino 
emissions up to about 100 TeV. As a matter of fact the neutrino interaction 
length is equal to the Earth diameter at about 40 TeV (Gandhi 2000). 
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The reader should notice that a rigorous calculation of event rates re- 
quires integrals since the neutrino flux, the effective area and the cross 
sections are functions of the energy. A useful parameter to calculate event 
rates is the effective area for i/'s, A e Jf (e.g. see Montaruli 2003, Ahrcns 
et al. 2004a), that is the sensitive area 'seen' by z/s producing detectable 
/i's when entering the Earth. In fact, the event rate in the detector from a 
point-like source at a declination 5 producing a differential flux of neutrinos 
is given by: 

= / dE v A e Jf(E v ,5)^-, (6) 

where E m i n and E max are the minimum and maximum energies of neutrinos 
for the considered flux, respectively. A e J* is given by: 

A e J f (E„, 8) = e - V gen ■ N AP <j v {E v ) ■ P E arth{E v ) (7) 

where e = jy"' (^P*^ . N gen is the number of generated events in the genera- 
tion volume Vg en (whose dimensions depend on the muon range) and N se i is 
the number of selected events that depends on track reconstruction quality 
cuts and selection criteria for background rejection; p ■ Na is the number of 
target nucleons per unit volume, with Na the Avogadro number and p the 
material density, a v is the neutrino CC cross section and 

P Earth (E u )^e- N ^yf^ dl (8) 

is the neutrino absorption probability in the Earth for a given neutrino of 
energy E v that propagates along paths of density p(l). The path in the 
Earth depends on the neutrino direction. 

The rejection of the residual atmospheric muons reaching detector depths 
can be achieved selecting upward-going muons induced by muon neutrinos 
interacting in the surrounding material. The rejection of backgrounds due 
to atmospheric neutrinos can be achieved through angular cuts around can- 
didate source directions, which mainly depend on the experimental angular 
resolution. For a telescope with sub-degree accuracy the number of atmo- 
spheric neutrinos is < f event /yr/km 2 . Also energy cuts can be applied 
to reduce further the backgrounds. Nevertheless, for galactic source models 
predicting neutrino emissions in the 1 — 100 TeV energy range, the tail of 
the atmospheric neutrino spectrum can still dominate over cosmic neutrino 
spectra. Hence angular cuts and possibly coincidence requirements with 
pulsed emissions or time modulations are more effective. 

Experimental results on neutrino searches from discrete galactic sources 
have been obtained by detectors of area < f 000 m 2 which have produced 
the first limits: 1MB (Svoboda et al. 1987; Becker-Szendy et al. 1995), 
Kamiokande (Oyama et al. 1989) and Kolar Gold Fields (Adarkar et al. 1991). 
Subsequently, limits were presented by the MACRO experiment (Ambrosio 
et al. 2001, updated in Montaruli 2003). Recently also Super-Kamiokande 
presented results for discrete steady sources (Washburn et al. 2003). 
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The new generation neutrino telescope AMANDA in the South Pole ice, 
with area larger than 10 4 m 2 , has not yet detected any excess from dis- 
crete sources, although it has set interesting limits for constraining models 
(Ahrcns et al. 2003; Ahrens et al. 2004a). AMANDA-II sensitivity, close 
to what expected for the ANTARES Mediterranean undersea neutrino tele- 
scope (Amram et al. 1999, Heijboer et al. 2003), could be sufficient to detect 
a signal among the atmospheric neutrino background for some of the models 
that are presented in this review (e.g. those in Sec.|2|and Sec. [BJ. Most of 
the models considered in this review predict neutrino event rates above the 
sensitivity limit of the future cubic-kilometer size detector IceCube (Ahrens 
et al. 2004b) and of a possible km 3 detector in the Mediterranean, as pro- 
posed to the European Community (KM3NeT 2004) and by the NEMO 
collaboration (Piattelli et al. 2003). 

Recent results on upper limits on the neutrino induced /i fluxes as a func- 
tion of declination for discrete sources are shown in Fig. (Montaruli 2003). 
In the region of declination observed by South Pole experiments, triangles 
indicate the upper limits obtained by AMANDA II (Ahrens et al. 2004a) 
using 699 upward events selected in 197 d during the year 2000. The most 
significant excess, observed around 21. lh R.A. and 68° declination, is of 8 
events and the expected background is 2.1. The probability to observe such 
an excess as a random fluctuation of the background is 51%. Three solid 
lines indicate the limits by AMANDA B-10 obtained in 130 days (Tab. 4 
in Ahrens et al. 2003), giving an idea of the effect of the statistical fluctua- 
tions of limits that are given in declination bands for different right ascen- 
sion bins. The lower line connects the lowest limits obtained in a declination 
band and the upper line connects the highest ones. The middle line connects 
the average of the limits in right ascension bins for each declination band. 
Limits from the experiments in the upper hemisphere are shown: Super- 
Kamiokande (Washburn et al. 2003, circles) and MACRO (Montaruli 2003, 
squares). Super-Kamiokande, with an angular resolution of about 2°, re- 
ported results using a sample of 2369 upward-going muons collected in 4.6 
yrs. The MACRO scintillator+tracking detector, with angular resolution 
< 1°, reported results using a sample of 1388 upward-going fi's collected 
in 6 yrs(Montaruli 2003). Ten neutrino events were measured inside a 3° 
half-width cone (that is expected to include 90% of an E~ 2 neutrino signal) 
around the plerion PSR1509-58, while 2 are expected from atmospheric i/'s. 
Even though PSR1509-58 is a source of interest as possible v emitter (see 
Sec. 0J, the significance of MACRO result is negligible when all the 1388 
directions of the measured events are looked at. Moreover, it is expected 
that E 1 "' 2 ^ 2 - 5 ) signals should produce at least 4-7 events in 1.5° around the 
source while only 1 is detected and the expected background is 0.5. A 7 
emission above 1.9 TeV was observed by CANGAROO from this source in 
1997 with 4cr significance, but not confirmed by the 1996 and 1998 data 
analysis with 2.5 TeV threshold (Sako et al. 2000). Models concerning this 
source are described in Sec.0]and predict much lower rates than what would 
be needed to explain the observed events in MACRO. From the same sources 
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Super-Kamiokande (A. Habig, private communication) observes the largest 
number of events between the selected catalogue looked at, but the data 
are still compatible with background fluctuations (9 events to be compared 
to a background of 5.4). The sensitivity of ANTARES (expected angular 
resolution of ~ 0.2° for E u > 10 TeV) is also shown for 1 yr of data taking 
(Heijboer et al. 2003). 




Fig. 1 90% c.l. upper limits on fi fluxes induced by v's with E 2 spectrum 
vs source declination for SK (green circles, Washburn et al. 2003), MACRO (red 
squares, Montaruli 2003, previous results in Ambrosio et al. 2001), AMANDA-B10 
(3 blue solid lines, Ahrens et al. 2003), AMANDA-II 2000 data (blue triangles, 
Ahrens et al. 2004a), ANTARES sensitivity after 1 yr (dotted black line, Heijboer 
et al. 2003). It has not been possible to apply a correction due to different [i average 
energy thresholds (SK ~ 3 GeV, MACRO ~ 1.5 GeV, AMANDA ~ 50 GeV). 
Nevertheless, the maximum of the response curves of these detectors for an E~ 2 
flux is at ~ 10 TeV, hence events contributing between 1-50 GeV should not 
make a large correction to these limits. 



3 Early phase of supernova explosion 

The supernovae can provide promising event rates during the relatively 
short period after the explosion. Nevertheless it should be considered that 
they are rare events, since the expected supernova rate in the Galaxy is of 
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the order of 2-4 per century and the pulsar birth rate is one per 20 -250 yrs 
(Lorimer 2003, Vranesevic et al. 2003). 

3.1 Supernova shock waves 

Type II supernovae are created by massive stars with extended envelopes. 
When the supernova shock breaks out of such a star, it creates a collision- 
less shock which may accelerate protons above 10 TeV (Waxman & Loeb 
2001). If the protons take a significant fraction of the post-shock energy, 
then observable fluxes of TeV neutrinos lasting ~ 1 hr can be produced in 
the inelastic collisions of relativistic protons with matter, about 10 hrs after 
the initial low energy neutrino flash following the core collapse of the star. 
The model predicts about 100 events in a 1 km 2 detector caused by muon 
neutrinos from a supernova at a distance of 10 kpc, if the proton acceleration 
efficiency is close to 10%. Clear predictions concerning the time lag between 
the low and high energy neutrino emission should allow to distinguish be- 
tween the proposed scenario and other models of neutrino production during 
the early phase of the supernova explosion. 

3.2 Supernovae with energetic pulsars 

Young supernova remnants (SNRs) have been proposed as sites of particle 
acceleration to high energy and, as a consequence, also possible sources 
of high-energy neutrino and 7-ray emissions. The first calculations of the 
neutrino flux from SNRs containing energetic pulsars have been performed 
by Sato (1977), Berezinsky & Prilutsky (1978), and more recently, also 
in relation to the explosion of SN 1981 A in the Large Magellanic Cloud 
(LMC), e.g. by Gaisser et al. (1987), Gaisser & Stanev (1987), Berezinsky 
& Ginzburg (1987), Berezinsky et al. (1988), Yamada et al. (1988), and 
Gaisser et al. (1989). 

In this review, we consider more recent models, starting from the one de- 
veloped by Protheroe et al. (1998). In this model the authors suppose that 
a pulsar with sufficiently short rotational period and high magnetic field 
has been formed just after the supernova explosion. Heavy ions, mainly 
56 Fe nuclei, can be accelerated from the star surface along open magnetic 
field lines into the pulsar magnetosphcrc. The slot gap model (see Arons 
& Scharlemann 1979) for the acceleration of nuclei has been applied. 5e Fe 
nuclei, extracted from the star surface by thermoionic emission, are gradu- 
ally accelerated in the slot gap and suffer photodisintegrations in collisions 
with the thermal photons coming from the neutron star surface. To check 
if photodisintegration of 56 Fe nuclei is important, Protheroe et al. (1998) 
compute the average number of neutrons extracted from 56 Fe nuclei as a 
function of their Lorentz factor, assuming that in each photodisintegration 
only one neutron is extracted. It is found that this process occurs efficiently 
only if the surface temperature is quite high, above « 10 7 K, i.e. within a 
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few days following the explosion. Some neutrons can be still extracted from 
the nuclei during the first year after the explosion, when the surface tem- 
perature drops to ~ 4 x 10 6 K (Nomoto & Tsuruta 1987). However, neutron 
star polar caps may have high temperatures, thus leading to efficient pho- 
todisintegration of iron nuclei emitted from the hot polar caps. The value 
of the temperature turns out to be dependent on the value of the magnetic 
field, the initial pulsar period and its radius. 

Relativistic neutrons extracted from Fe nuclei in this way move bal- 
listically through the magnctosphcrc and beyond the light cylinder. Their 
Lorentz factors, which lie within the range ~ 10 4 - 10 7 , depend on the 
pulsar parameters and on the temperature of the radiation. Typical values 
of the temperature are of the order of ~ 10 6 — 10 7 K for periods shorter 
than ~ 10 ms and surface magnetic fields 10 12 G. Such neutrons decay at 
a distance of ~ 0.1 — 100 pc from the pulsar. Initially, when the SNR shell 
is still opaque to neutrons, they mostly interact with target nuclei as they 
travel out through the dense SNR shell, thus producing 7-rays and neutrino 
signals. 

Protheroe et al. calculate the (y^ + v^) neutrino spectra for a source 
located at a distance d — 10 kpc, at a time t = 0.1 year after the explosion, 
assuming initial periods of the pulsar Pq = 5 ms and Pq = 10 ms. The mag- 
netic field is kept fixed at B = 10 12 G. Since the ion emission is dependent 
on the temperature of the star surface and the polar cap, two limiting cases 
have been considered, i.e. 

i) " no polar cap heating" , where the whole star has temperature equal to 
the surface temperature; 

ii) "maximum polar cap heating", where the polar cap has temperature 
larger than the surface temperature. 

In the following, we report the neutrino energy spectra, calculated in units 
of GeV -1 m~ 2 s _1 , for an initial pulsar period Po = 10 ms (see Fig. 9 in 
Protheroe et al. 1998) 





3.622 



(E„ > 10 4 GeV) (no h.) 



(9) 




(10) 
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For an initial period P = 5 ms, the neutrino spectra are: 



-1.480 



/ JP \ -3-687 

= 1.126 x 10 6 ( J (E„ > 10 4 GeV) (no h.) 

(11) 



/ TP \ - 1 ' 788 

9.120 x 10- 4 f YGeV ) ^ < 1()4 GeV ) 

(jp \ -3.687 
j^-j J (E„ > 10 4 GeV) (max. h.) 

(12) 

From these calculations it is clear that selecting an angular region of 1° 
around the source, the neutrino flux above 1 TeV is well above the atmo- 
spheric neutrino background. For a 10° region, a neutrino signal from a 
galactic source for the Po = 5 ms case may be observable above the back- 
ground at energies above 10 TeV. However, as shown in Protheroe et al. 
(1998), the neutrino light curve drops rapidly at t > 0.1 year, and the sig- 
nal could be detected among the atmospheric background only in detectors 
of dimensions of the order of km 2 . 

Since the neutrino emission from acceleration and subsequent fragmen- 
tation of nuclei is beamed, the neutrino signal should be strongly anisotropic 
and modulated with the pulsar period. Hence the signal could be observed 
only if the beam is in the Earth direction. 

Another target for hadrons accelerated by a pulsar early after the super- 
nova explosion can be provided by the thermal emission of the expanding 
opaque supernova envelope. Beall & Bednarek (2002) consider neutrino pro- 
duction in collisions of nuclei with the radiation filling the cavity below the 
hot expanding supernova remnant envelope. The heavy nuclei, accelerated 
close to the light cylinder of the pulsar or in the pulsar wind region, are 
completely disintegrated in collisions with the dense radiation field at early 
times after the supernova explosion. The resulting products (protons and 
neutrons) lose energy efficiently in pion production. However, due to the 
very dense radiation field, pions can also lose energy before decaying due 
to Inverse Compton Scattering (ICS). These ICS energy losses of pions are 
important during the first ~ 10 4 s after supernova explosion. However, later 
pions decay into muon neutrinos which escape without absorption through 
the expanding supernova envelope. When the column density of matter 
through the envelope drops to ~ 10 3 g cm -3 , the radiation is not further 
confined by the envelope and the nuclei interact with the envelope matter 
producing muon neutrinos. Muon neutrino fluxes have been calculated from 
interactions of nuclei with the radiation field inside the envelope during the 
period 10 4 -j- 2 • 10 6 s after the supernova explosion and from interactions 



dN v 
dE v 
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with the matter envelope during the period 2 • 10 6 -=-3 • 10 7 s after the explo- 
sion. During these periods the density is low enough so that pions produced 
in the interactions can decay. At later times the neutrino flux drops signifi- 
cantly because the column density of matter through the envelope drops to 
several g cm -3 and particles are not completely cooled in collisions, hence 
they do not lose energy efficiently. 

The model predicts up to ~ 10 4 — 10 5 events in a 1 km 2 neutrino detector 
from a supernova at a distance of 10 kpc during the first year after the 
explosion, for pulsars with initial period of the order of a few milliseconds 
and typical surface magnetic fields of a few 10 12 G, or pulsars rotating with 
a period of 10 ms and with surface magnetic field of the order of 10 14 G 
provided that the acceleration efficiency of nuclei by the pulsar is £ = 1, i.e. 
all rotational energy of the pulsar goes into nuclei acceleration. Hence, these 
results hold if these pulsars are able to accelerate nuclei up to energies of the 
order of 10 20 eV. The event rates, obtained by using the neutrino detection 
probability in Gaisser & Grillo (1987) and the absorption coefficients of 
neutrinos in the Earth calculated by Gandhi (2000), are given in Tab. 0] 
of Beall and Bednarek (2002) for £ = 1. The neutrinos are produced with 
a E~ x power law differential spectrum between ~ 0.1 — 100 PeV for the 
case of hadron interactions with radiation, and up to ~ 10 PeV for hadron 
interactions with the envelope matter (see Figs. 2 and 3 in Beall & Bednarek 
2002). 

However, the value of the acceleration efficiency £ is a free parameter 
of the model, and might be constrained by observations. For this purpose, 
Blasi et al. (2000) propose that heavy nuclei injected by fast pulsars into the 
Galaxy, after the supernova becomes transparent, can explain the observed 
cosmic ray (CR) flux at ~ 10 20 eV. In fact, from the normalization of the 
flux of nuclei predicted in this model to the observed cosmic ray flux in the 
Galaxy (see Eq. 8 in Blasi et al. 2000), the acceleration efficiency of a pulsar 
can be estimated, and the following relation can be obtained 

£eQ/(r 2 i? 2 B 13 ) «4x HT 6 , (13) 

where e is the fraction of pulsars able to accelerate iron nuclei to 10 20 eV, 
Q is the trapping factor of nuclei in the Galactic halo, r = 100t2 is the rate 
of neutron star formation, R = 10i?i kpc is the radius of the Galactic halo, 
and B = 10 13 £?i3 G is the surface magnetic field of the pulsar. Applying 
reasonable values to these parameters, Q = 1, r 2 = 1, R\ = 3, e = 0.1, 
and Bi 3 also considered in the pulsar model, the estimated efficiencies of 
pulsar acceleration £ are 3.6 • 10 -3 for B%2 = 100 and 1.4 • 10 -4 for B%2 = 4. 
After taking into account these values of £ the event rates are shown in 
the first two columns of Tab. 0] for the indicated pulsar parameters. The 
muon neutrino event rates expected for horizontal and vertical directions, 
after including the absorption by the Earth, are indicated by (H) and (V), 
respectively. The reactions N-7 — > and Fe-M^ are used for the process 
of neutrino production in hadron (nucleons or iron nuclei) collisions with the 
radiation and with the envelope matter, respectively. The model predicts up 
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Table 1 (y^ + induced events from a supernova at distance of 10 kpc in the 
horizontal (H) and vertical (V) directions in a km 2 detector. The calculation is 
done for different pulsar periods and magnetic fields (from Tab. 1 in Beall and 
Bednarek (2002)). The calculation of rates includes the indicated values of the 
acceleration efficiency of the pulsar £ (see eq. 1130 . 



pulsar =>■ 


-Puis — 3 


P ms = 10 


P ms = 20 


process Jj- 


B 12 =4 


Bia = 100 


B 12 =4 


N-7 - (H) 


1.2 


86.4 


2.6 


Fe-M^ u M (H) 


12.2 


648.0 


11.3 


N-7 -» (V) 


0.3 


18.4 


1.1 


Fe-M-> v M (V) 


4.8 


230.4 


8.6 


t 


1.4 x 10" 4 


3.6 x 10" 3 


1 



to hundreds of neutrino events in a 1 km 2 detector for pulsars, with initial 
periods of several ms and surface magnetic fields of the order of a 10 14 G,at 
a distance of 10 kpc provided that £ > 10 -3 . On the other hand, Crab- 
like pulsars at a distance of 10 kpc, with initial periods of several ms and 
surface magnetic fields of the order of a few 10 12 G, can produce detectable 
event rates in a 1 km 2 detector provided that £ is not far from 1 (see third 
column in Tab. 0. As a matter of fact, a pulsar with a period of 20 ms 
and a magnetic field of 4 ■ 10 12 G cannot accelerate nuclei up to 10 20 eV, 
hence the acceleration efficiency is not constrained by eq. 1131 Hence, as an 
example, event rates are given in the optimistic assumption of fully efficient 
pulsar acceleration. 

Neutrino production in high energy hadron-hadron collisions in young 
pulsar driven supernovae has been further investigated by Nagataki (2004). 
It is assumed that most of the rotational energy of a pulsar with millisecond 
initial period is taken by relativistic hadrons with the Lorentz factor kept 
as a free parameter of the model. Hadrons, accelerated in the pulsar wind 
region, are thermalized in the inner shock and interact inside the pulsar wind 
nebula between themselves. The neutrino event rates with energies above 
10 GeV expected from a 1 ms pulsar at a distance of 10 kpc, are above the 
atmospheric neutrino background, and could be detected by cubic-kilometer 
neutrino telescopes during the first 100 yrs after the explosion. However, 
note that according to this modeling the shock in the pulsar wind is located 
very close to the pulsar, at a distance of ~ 10 12 cm (for a pulsar with 
initial period 1 ms, surface magnetic field 10 12 G, and age of 100 yrs). 
Hence, the density of relativistic hadrons can be high enough for efficient 
production of pions, which subsequently decay into 7-rays and neutrinos. 
Nevertheless, it seems that this model may have difficulties in predicting the 
observed location of the shock in the Crab Nebula, which is at a distance 
of ~ 3 x 10 17 cm at its present age of ~ 10 3 yrs. Unless this constraint is 
fulfilled, the predicted fluxes of neutrinos could be not reliable, since, as the 
author shows in the paper, the expected neutrino flux is very sensitive to 
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the shock location. Moreover, the amount of energy supplied by the pulsar 
to the expanding nebula is not taken into account. For a pulsar with 1 ms 
initial period and magnetic field 10 12 G, the energy of relativistic protons, 
which is equal to the rotational energy lost by the pulsar during the first 10 
yrs (~ 10 52 ergs) , is much larger than the initial kinetic energy of a typical 
supernova envelope (~ 10 50 ergs for the Crab Nebula). This additional 
source of energy should accelerate the nebula, moving the location of the 
shock outwards. As a result, the relativistic hadron densities decrease and 
the expected neutrino fluxes should be much smaller. 

4 Pulsar wind nebulae (Plerions) 

A number of plerions have been discovered in radio, optical, and X-ray 
bands, with the Crab as the youngest and most energetic source. Only 
about 10% of supernova remnants are classified as plerions. Plcrionic na- 
ture is usually indicated by a center-filled morphology, resulting from the 
continuous injection of pulsar electrons into the nebula, with the additional 
constraint that the spectra must be non-thermal (but power-laws) resulting 
from a statistical acceleration process. The radio, optical and X-ray emis- 
sion observed from plerions is believed to be due to synchrotron emission. 
ICS is usually seen at higher energies compared to synchrotron emission. It 
is likely that hadrons can also play an important role in radiation processes 
(e.g. Cheng et al. 1990, Aharonian & Atoyan 1995). 

In order to be more quantitative, let us suppose that TeV photons pro- 
duced by a plerion are the decay products of 7r°'s. Neutral pions can be 
produced by accelerated protons in either pp or p-f collisions, depending on 
the relative target density of photons and protons in the source region where 
protons are accelerated. The {v^ + v^) neutrino flux, dN„/dE„, produced by 
the decay of charged pions in the source can be derived from the observed 7- 
ray flux, dN 7 /dE 7 , by imposing energy conservation and neglecting photon 
absorption: 



where E™ m (E™ ax ) is the minimum (maximum) energy of the photons that 
have hadronic origin. E™ m and E™ ax are the corresponding minimum and 
maximum energies of the neutrinos. The factor K depends on whether the 
7r°'s are of pp or p-f origin. As a first order approximation, using average 
values, it is generally assumed that K = 1 for pp scattering, whereas K = 4 
for P7 interactions. 

Such simple rescaling of the 7-ray fluxes to the neutrino fluxes has been 
done by Guetta and Amato (2003) for a few pulsar wind nebulae (PWNe) 
observed in TeV 7-rays, i.e. the Crab Nebula, the Vela SNR, the pulsar 
wind nebula around PSR1706-44 and the nebula surrounding PSR1509-58. 
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For instance, the neutrino flux predicted for the Crab nebula (Guetta 2003, 
private communication) is given (in units of GeV -1 m~ 2 s _1 ) by: 



It is concluded that all these PWNe could be detected by a 1 km 2 neutrino 
detector with event rates per year between 1 up to about 12 in the energy 
range E v ~ 1 — 100 TeV. However, these estimations are based on the 
assumption that all 7-rays above 2 TeV are produced in pp collisions via 
pion decay. It has been argued by Bednarek & Bartosik (2003) that this is 
a reasonable approximation for relatively young nebulae, of the age of the 
Crab nebula or younger, but not for older nebulae, such as those mentioned 
above. 

A more specific model for 7-ray and neutrino production in the Crab 
Nebula has been proposed a few years ago by Bednarek and Protheroe 
(1997). The general scenario is similar to the one previously discussed for 
the very young supernova remnants (Protheroe et al. 1998): iron nuclei from 
the Crab pulsar are accelerated in the outer gap (Cheng et al. 1986) and 
photodisintegrate in collisions with the non-thermal radiation field filling 
the outer gap, thus producing energetic neutrons which are injected into the 
nebula. Neutrons can decay either inside or outside the Crab nebula. Those 
neutrons which decay inside the nebula produce protons, which collide with 
the matter in the nebula producing 7-rays and neutrinos. Neutrons which 
decay outside the nebula can be captured by the expanding nebula at later 
times or contribute to galactic cosmic rays. The neutrino fluxes predicted by 
this model are dependent on the Crab pulsar and nebula parameters such as: 
the initial pulsar period, the nebula radius and the mass ejected during the 
supernova explosion. When reasonable values are taken into calculations, 
then the expected neutrino spectrum (in units of GeV -1 m~ 2 s _1 ) can be 
described by the following broken power law: 



If the observed TeV 7-ray emission from the Crab Nebula is due to the 
hadronic processes at energies above 10 TeV, then the model predicts a few 
neutrino events per year in a 1 km 2 detector. 

Recently, two papers discuss independently the hadronic processes in the 
PWNe, assuming that a significant part of the rotational energy lost by the 
pulsar is taken by relativistic heavy nuclei accelerated in the pulsar wind 
region (Bednarek & Bartosik 2003, Amato et al. 2003). This assumption 
can explain some fine structures in the Crab Nebula, as discussed by, e.g., 
Gallant & Arons (1994). According to their model, the nuclei can generate 
Alfven waves just above the pulsar wind shock, which resonantly scatter 
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off positrons and electrons, accelerating them to the energies observed, e.g., 
in the Crab Nebula. Relativistic nuclei injected into the nebula can inter- 
act with the nebula matter, and produce 7-rays and neutrinos via pion 
decay. The consequence of this hybrid hadronic-leptonic model is that it 
self-consistently treats both leptonic and hadronic processes. 

In the model proposed by Bednarek & Bartosik (2003), the basic param- 
eters of the expanding nebula, the injection spectra of hadrons and leptons, 
and their equilibrium spectra inside the nebula, depend on time, according 
to the expectations of the pulsar evolution model. Some free parameters of 
the model are fixed by the comparison with the TeV 7-ray emissions from 
PWNe, such as the Crab Nebula, the Vela nebula, and the nebula around 
PSR 1706-44. Based on that, they predict the TeV 7-ray fluxes from PWNe 
which contain energetic pulsars. The correct modeling of the known X-ray 
emission from the TeV 7-ray PWNe is possible if it is assumed that the 
lepton acceleration efficiency by nuclei drops with the age of the nebula. 
This suggests that the lepton acceleration efficiency for a specific nebula 
can also depend on time. This effect has not been included into the calcula- 
tions. In this model, Bednarek (2003a) calculates the neutrino spectra from 
a few PWNe and predicts also the neutrino event rates in a 1 km 2 detector 
from a few PWNe. It is concluded that the largest rate is from the Crab 
Nebula that would produce a rate larger than 1 event per year (using the 
detection probabilities by Gaisser & Grillo 1987). In this case, the predicted 
spectrum is above the atmospheric neutrino background in the energy range 
- 10 - 100 TeV, and is given (in units of GeV" 1 m~ 2 s" 1 ) by: 



Other PWNe can be detected only if they are close to high density regions, 
which is probably the case of MSH15-52 around PSR1509-58. The ratio of 
the TeV neutrino and 7-ray fluxes depends on the age of the PWNa in 
this hybrid hadronic-leptonic model, dropping from ~ 58% for the Crab 
Nebula, down to 22% for the Vela nebula, and 0.5% for the nebula around 
PSR 1706-44. 

Amato et al. (2003) present an analytical approach to the leptonic and 
hadronic processes inside the PWNe. The neutrino event rates predicted 
by them for the Crab Nebula in a 1 km 2 neutrino detector are between 4 
and 14, depending on the Lorentz factor of nuclei injected by the pulsar. 
They are significantly higher than those obtained by Bednarek (2003a). The 
discrepancies between the models might be ascribed, on one hand, to the 
different treatment of the injection spectra of hadrons, which actually evolve 
in time. On the other hand, the neutrino spectra are determined by the 
hadronic equilibrium spectra inside the nebula, and these are constrained 
by the observations. A time dependent treatment of the equilibrium spectra 
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is required in order to have the correct neutrino spectrum. In the model by 
Amato et al. (2003), the Lorentz factors of the nuclei are assumed to be 
constant and independent on the pulsar parameters, which actually should 
evolve in time due to the pulsar rotational energy losses. 

5 Shell-type supernova remnants 

A few spherical supernova remnants, such as SN1006 (Tanimori et al. 1998), 
Cas A (Aharonian et al. 2001), and SNR RX J1713. 7-3946 (Enomoto et 
al. 2002), have been reported to emit 7-rays with energies above ~ 1 TeV. 
From these observations, it cannot be established whether the mechanism 
of 7-ray production in these sources is leptonic or hadronic. In other words, 
these high-energy photons could be produced either by electron beams via 
ICS (Pohl 2001), or by proton beams via the production and subsequent 
decay of neutral pions (Gaisser et al. 1995; Pohl 2001). In the latter case, the 
high-energy photons are accompanied by neutrinos because charged pions 
are produced along with 7r°'s. 

Recently the CANGAROO collaboration has shown that the energy 
spectrum of the 7-ray emission from the supernova remnant RX J1713.7- 
3946 matches that expected if the 7-ray are the products of 7r° decay gener- 
ated in pp collisions (Enomoto et al. 2002). Alvarez-Muniz & Halzen (2002) 
have estimated the expected rate of neutrino events in a 1 km 2 detector 
using the normalization between 7-ray and neutrino fluxes (see eq. I14fl . 
Assuming that all charged pions decay in the environment of the super- 
nova and that the neutrino spectrum follows the input proton spectrum, 
assumed to be of ~ E~ 2 type, the calculated (y^ + v^) flux is given (in 
units of GeV -1 m~ 2 s _1 ) by: 



Alvarez-Muniz & Halzen have predicted large muon neutrino event rates 
from SNR RX J1713. 7-3946, of the order of - 40 km" 2 yr" 1 . However, the 
reader should notice that the hadronic origin of 7-rays from this supernova 
remnant has not been confirmed by the analysis of hadronic and leptonic 
processes performed by Reimer & Pohl (2002). 

Neutrino fluxes have been predicted by Alvarez-Muniz & Halzen for 
some other supernovae, such as, Cassiopeia A, Crab Nebula, and the shell- 
type supernova remnant Sagittarius A East at the Galactic Centre (Buckley 
et al. 1997). For the SNR Cassiopeia A, the predicted neutrino event rate 
is ~ 5 km -2 yr _1 , whereas for Sagittarius A East it is much larger (~ 
140 km -2 yr _1 ). As far as the Crab nebula is concerning, the event rates 
turn out to be close to the ones expected from the atmospheric neutrino 
background. With these predicted fluxes, RX J1713. 7-3946 and Sagittarius 
A East are sources of interest also for telescopes of smaller dimensions than 1 
km 2 , such as the ANTARES telescope (Amram et al. 1999) and AMANDA- 
II (Ahrens et al. 2004a). 
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6 Pulsars in high density regions 

The lifetime of massive stars which explode as supernovae, with subsequent 
formation of pulsars, is relatively short, of the order of 10 6 — 10 7 yrs. There- 
fore, pulsars are likely to appear inside or close to massive stellar associations 
where the amount of distributed matter is very high. Particles accelerated by 
pulsars can be trapped by this high density medium, which is characterized 
by magnetic fields of the order of 10~ 5 — 10~ 4 G, or even ~ 1CP 3 G in- 
side huge molecular clouds (Crutcher 1999). The largest stellar associations 
in our Galaxy contain hundreds early type massive (OB) and Wolf-Rayet 
(WR) stars. Hence, the average formation rate of pulsars inside the massive 
associations can be estimated to be about 10 -4 per year. Therefore, it is 
likely that on average a relatively young pulsar, with an age of ~ 10 4 yrs, 
should be present in every massive association. Such pulsars should be ef- 
ficient accelerators of particles up to very high energies, as it is proved by 
their pulsed 7-ray emission and surrounding non-thermal nebulae. 

6.1 The Galactic Centre dense region 

The most interesting site for the occurrence of the processes described above 
is the huge concentration of mass in the inner ~ 50 pc around the Galactic 
Centre (GC), with total mass of ~ 10 6 M Q , average density 10 2 cm~ 3 , and 
magnetic field 3 x 10~ 5 G. This region contains also a huge molecular cloud 
with radius 10 pc, density 10 3 cm~ 3 , and magnetic field 10~ 4 G. The GC 
is also a site where high energy processes occur. For example, the EGRET 
detector on board of the Compton Gamma-Ray Observatory satellite has 
detected an intense 7-ray source with luminosity 2 x 10 37 erg s -1 (Mayer- 
Hasselwander et al. 1998). There are also some indications of TeV 7-ray 
emissions from this region (Kosack et al. 2003, Tsuchiya et al. 2003). The 
AGAS A collaboration reported the existence of an extended excess of cosmic 
rays (CR) of energy around ~ 10 18 eV from the direction close to the GC 
(Hayashida et al. 1999). This excess has been confirmed by the analysis of 
the SUGAR data (Bellido et al. 2001). 

Recently, Bednarek (2002) has considered a scenario in which energetic 
pulsars accelerate nuclei inside (or close to) this high density region to en- 
ergies observed by the above mentioned cosmic ray experiments. Such a 
pulsar should have an initial period P = 10~ 3 -Po, ms s and surface magnetic 
field B = 10 12 i?i2 G, which fulfill the condition Bi 2 Po, ms > 0.4 in order to 
accelerate hadrons above 10 18 eV. As an example, a pulsar with P = 2 ms 
and B — 6 x 10 12 G has been chosen. The nuclei accelerated by the pulsar, 
as discussed in the model by Bednarek & Bartosik (2003), escape from the 
pulsar nebula and are partially captured by the high density region in the 
GC. They partially disintegrate in collisions with matter, thus injecting high 
energy neutrons which are responsible for the observed excess of cosmic rays 
with energy ~ 10 18 eV. Moreover, these nuclei produce 7-rays and neutri- 
nos. Bednarek (2002) has also investigated the likely parameter range for 
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Table 2 event rate per year in a km 2 detector from the Galactic Centre for 
an E' 2 - 6 spectrum between 10-100 TeV and ~ 10 PeV 



medium =>■ 
pulsar age J| 


GC extended 
region 


dense molec. 
cloud 


3 x 10 3 yrs (H - V) 


3.8 - 2.0 


16 - 11 


10 4 yrs (H - V) 


8.8 - 5.3 


30 - 23 



the GC medium, considering the case of hadrons trapped inside the dense 
molecular cloud with density 10 3 cm -3 , radius 10 pc, and magnetic field 
10 -4 G, and a larger region around the GC with density 10 2 cm -3 , radius 
50 pc, and magnetic field 3 x 10~ 5 G. From the normalization of the calcu- 
lated neutron flux to the observed flux of CR particles in excess, he predicts 
the fluxes of TeV 7-rays, which should be detected by the next generation 
of Cherenkov experiments, and the fluxes of neutrinos at the Earth. The 
neutrino spectra are characterized by power laws with a spectral index of 
~ —2.5 from ~ 10 — 100 TeV up to ~ 10 PeV. It is found, considering 
the probability of detecting muons produced by such neutrinos (Gaisser & 
Grillo 1987) and the neutrino absorption by the Earth (Gandhi 2000), that 
a few up to several muon neutrino events should be observed in a 1 km 2 
neutrino detector per year for the horizontal (H) and the vertical (V) direc- 
tions from a pulsar which was formed (3 — 10) x 10 3 yrs ago (see Tab. |3 for 
the details). 

The observed CR excess at ~ 10 18 eV is also interpreted as due to 
the last (or last few) Gamma Ray Bursts (GRBs) which occurred in our 
Galaxy (Biermann et al. 2004). If the rate of GRBs is of the order of one 
per 10 6 yrs, then hadrons accelerated by the GRB relativistic shock convert 
partially to neutrons which decay into protons inside the Galaxy, and after 
some time protons are partially converted into neutrons in collisions with 
matter. The relative flux of neutrinos to the flux of neutrons, produced in 
these hadronic collisions, with energies below ~ 10 18 eV has been estimated. 
It should produce similar event rates in a 1 km 2 neutrino detector as those 
estimated by Bednarek (2002). 

The relativistic hadrons can also be accelerated by the shock waves of 
the supernovae which appeared in the central high density region. One of 
such shell-type SNR associated with Sgr A East, with an age of ~ 10 4 yrs, 
has been recently considered as responsible for the unidentified EGRET 
source at the GC (see Fatuzzo & Melia 2003). These relativistic hadrons 
should interact with the high density medium (10 3 cm -3 ) producing 7-rays 
and neutrinos. However, expected neutrino fluxes at TeV energies have not 
been estimated yet in the context of this model. 
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6.2 The massive stellar associations 

The largest massive stellar association, at a relatively small distance of 1.7 
kpc, is Cygnus OB2. Its total mass derived from the CO survey is ~ 3.3 x 10 5 
M Q (Butt et al. 2003) and the total stellar mass is (4 - 10) x 10 4 M Q 
(Knodlsedler 2000). Moreover, it contains many massive stars of the O type, 
- (120 ± 20) (Knodlsedler 2000). Recently, Cygnus OB2 became also an 
interesting high energy source due to the reports of detection of a few 7-ray 
EGRET sources (e.g. Mori et al. 1997, Lamb & Macomb 1997), a TeV 7-ray 
source by HEGRA (Aharonian et al. 2002), and an excess of CR particles 
at energies ~ 10 18 eV by AGASA (Hayashida et al. 1999). Therefore, it is 
also an interesting potential source of high energy neutrinos. 

Torres et al. (2004) propose that TeV 7-rays and neutrinos can be pro- 
duced by CRs through interactions with the nuclei in the innermost parts 
of the winds of massive O and B stars. This model would explain the TeV 
emission observed by HEGRA in the vicinity of Cygnus OB2 (Aharonian ct 
al. 2002). The main acceleration region of CRs would be the outer boundary 
of the supper-bubble produced by the core of the association of OB stars. 
Moreover, only CRs with energies > 1 TeV are able to penetrate dense inner 
regions of the winds from outside since for them the diffusion timescale in 
the wind is shorter than the convection time scale for likely parameters of 
the massive stars. This fact explains why a GeV 7-ray source related to the 
reported TeV 7-ray source by HEGRA has not yet been observed. On the 
other hand, a nearby EGRET source, 3EG J2033+4118, is coincident with 
the center of the association, where the GeV emission might be produced in 
the terminal shocks of powerful stars existing therein, in the colliding wind 
binary system Cyg OB2 No. 5 or in a combination of these scenarios. In or- 
der to explain the observed flux from the unidentified TeV 7-ray source, the 
model requires an enhancement of the CR flux in the region of Cygnus OB2 
of the order of 300 compared to the CR flux in the vicinity of the Earth, if 
the ISM density is of the order of 0.1 cm -3 . The predicted neutrino signal 
from Cyg OB2 above 1 TeV in the IceCube neutrino detector should be 
a factor ~ 1.8 larger than the atmospheric neutrino background in search 
bins of 1° x 1°. 

Another model able to explain possible high energy phenomena in Cygnus 
OB2 assumes the birth of a very energetic pulsar about 10 4 yrs ago inside 
this association (Bednarek 2003b) . Taking into account the number of mas- 
sive stars inside Cygnus OB2 (~ 100) and their characteristic lifetime (of a 
few 10 6 yrs), such pulsars are likely to be produced with the required rate 
of one per ~ 10 4 yrs. This pulsar, with reasonable parameters, such as an 
initial period of 2 ms and a surface magnetic field of 6 x 10 12 G, would 
be able to accelerate heavy nuclei to Lorentz factors of ~ 10 9 soon after 
formation. Part of these nuclei can be captured for a long time in the dense 
regions of Cygnus OB2, characterized by a strong magnetic field (~ 10~ 4 
G), producing relativistic neutrons in collisions with matter. These neu- 
trons can be responsible for the observed CR excess observed in AGASA 
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(Hayashida et al. 1999). On the other hand, the pulsed emission of this Vela 
type pulsar, ~ 10 4 yrs old, could be responsible for the observed GeV 7-ray 
source. In fact, the level of 7-ray emission and the shape of the spectrum of 
the Vela and PSR 1706-44 pulsars are consistent with the observations of 
the EGRET source 3EG J2033+4118 in the direction of Cyg OB2. The TeV 
7-ray emission from the nebula surrounding this pulsar can in turn explain 
the HEGRA TeV 7-ray source (see for detailed modeling Bednarek & Bar- 
tosik 2003). However, the model considered by Bednarek (2003b) predicts 
that the fluxes of muon neutrinos, produced by nuclei within the pulsar 
nebula which in turn are captured by the massive association, are too low 
to be detected by a 1 km 2 neutrino detector (~ 0.5 neutrino events from 
the pulsar nebula and ~ 0.14 from the dense region of the Cygnus OB2 per 
year). 

It has also been suggested that rclativistic nuclei may suffer strong 
photo-disintegrations in the radiation field present inside the massive as- 
sociation Cyg OB2 (Anchordoqui et al. 2003a). This soft photons are pro- 
duced by the above mentioned population of the massive stars and the warm 
molecular clouds. Neutrons extracted from the nuclei, with energies below 
~ 2 x 10 eV, decay into antineutrinos during their propagation towards 
the Earth along distances smaller than 1.7 kpc. A few antineutrinos of all 
flavors (since the electron antineutrinos from neutron decay are expected to 
oscillate into other flavors) might be detectable per year by a 1 km 2 detector 
if the spectrum of nuclei inside the Cyg OB2 is very steep oc E~ 3A in the 
PeV-EeV region (Anchordoqui et al. 2003a). Nevertheless these calculations 
do not consider that the interactions of nuclei with matter (as considered by 
Bednarek 2003b) might be more important than those with radiation. If the 
total amount of gas inside the core of Cyg OB2 (within 10 pc) is <~ 3 x 10 4 
M Q (corresponding to an average density of ~ 300 cm~ 3 ), then neutrinos 
produced in collisions of nuclei with matter would dominate over neutrinos 
from decays of neutrons, produced in photodisintcgrations of nuclei. This 
would be due to the larger (by about an order of magnitude) cross sections 
for collisions of nuclei with matter, with respect to radiation. 



7 Neutron stars in binary systems 

The calculations of neutrino fluxes from binary systems have been reviewed 
by Gaisser et al. (1995). Most of the earlier models assumed that particles 
are accelerated due to the energy generated during the accretion of matter 
onto a slowly rotating neutron star (e.g. Gaisser & Stanev 1985, Kolb et 
al. 1985, Berezinsky et al. 1985). Those models predict high neutrino fluxes 
since they are normalized to the measurements of VHE 7-ray emissions from 
Cyg X-3, still nowadays not well understood. In another model, Harding & 
Gaisser (1990) consider hadrons accelerated in the pulsar wind shock taking 
as a reference the massive binary Cyg X-3. This model predicts detectable 
fluxes of neutrinos in a 0.1 km 2 detector (Gaisser et al. 1995) only for the 
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case of eclipsing binaries, when the massive star is between the source of 
relativistic hadrons and the observer. 

A similar scenario has been more recently considered by Bartosik et 
al. (2003), i.e. a binary system containing an energetic pulsar and a mas- 
sive stellar companion, with the purpose to calculate the neutrino fluxes 
produced by pulsar accelerated nuclei. Provided that the binary system is 
compact enough and the star is luminous, the nuclei suffer photodisinte- 
grations in the thermal radiation field of the massive star. The neutrons 
from disintegration move ballistically, some of them interact with the mas- 
sive star matter. On the contrary, the nuclei which survived disintegration 
and the secondary protons are injected through the shock structure into the 
magnetic field in the massive star wind region. Bartosik et al. calculate the 
paths of charged hadrons in the massive star magnetic field. It is found that 
some of these charged hadrons can fall at large angles (larger than expected 
from the simple linear propagation) onto the massive star, thus produc- 
ing neutrinos. The neutrino fluxes and spectra produced by neutrons and 
charged hadrons are calculated for different locations of the observer with 
respect to the orbital plane of the binary system. It is found that interest- 
ing fluxes of neutrinos are emitted at angles larger than those intercepted 
by the massive star. Therefore, not only eclipsing binaries might become 
detectable sources of neutrinos by the future telescopes. 

Slowly rotating, accreting neutron stars have been also considered re- 
cently as a source of neutrinos. It is proposed that if a neutron star is 
strongly magnetized and spins down slower than the accretion disk, then 
an accelerating gap should appear in the pulsar magnetosphere above the 
accretion disk (Cheng & Ruderman 1989). Hadrons can be accelerated in 
such a gap up to ~ 100 TeV, if the disk luminosity is of the order of 10 37 
erg s _1 , and the surface magnetic field of the neutron star close to 10 12 
G. In the frame of this model, Anchordoqui et al. (2003b) estimate the 
neutrino flux produced in cascades initiated by relativistic hadrons in the 
matter of the optically thick accretion disk. The authors predict a neutrino 
signal above the atmospheric neutrino background in a region of 1° x 1° 
around the source (signal to noise ratio w 1.7) between 300 GeV and 1 TeV 
from the binary system A0535+26. The predicted neutrino flux, in units of 
GeV -1 rar 2 s _1 , is given by 



The signal should appear when the neutron star is close to the periastron 
due to the formation of a dense accretion disk. Therefore, it should be 
modulated with the period of the binary system. The neutrino signal from 
A0535+26 should be above the sensitivity of the IceCube detector during a 
few binary system reversals. 
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8 Microquasars 

Microquasars are galactic X-ray binary systems, which exhibit relativistic 
radio jets. Therefore, similar processes of neutrino production by relativistic 
particles as occurring in active galactic nuclei (AGNs) can be considered, 
e.g. proton initiated cascade model in the jet (Mannheim & Biermann 1992, 
Miicke & Protheroe 2001), in the radiation of an accretion disk (Bednarek & 
Protheroe 1999), or interaction of neutrons with the matter of an accretion 
disk (Nellen et al. 1993). 

Recently, Levinson & Waxman (2001) have proposed a model for neu- 
trino production in microquasars, similar to the AGN proton initiated cas- 
cade model of Mannheim et a. (1992). The authors argue that hadrons can 
be accelerated up to ~ 10 16 eV in the inner part of the jet by internal shocks 
in the jet. Pions are produced in collisions of hadrons with the external X- 
ray photons, and with the synchrotron photons from the jet produced by 
leptons accelerated in this same shock. The neutrino fluxes, expected in the 
energy range 1-100 TcV, are very high, especially from persistent sources, 
e.g., ~ 10 3 neutrino event rates in a 1 km 2 neutrino detector in the case 
of SS433. For the case of bursting sources, such as GRS 1915+105, the ex- 
pected sensitivity of a 1 km 2 detector would require observation of about 30 
outbursts for detection. In a subsequent paper, Distefano et al. (2002) esti- 
mated the expected neutrino event rates in a 1 km 2 detector from the large 
population of known microquasars, assuming that relativistic protons take 
10% of the jet power. The authors conclude that persistent microquasars, 
such as SS433 or GX 339-4, should produce more than 100 neutrino events 
in such a 1 km 2 detector. Therefore, some of the neutrino events could 
be observed by AMANDA-II. Also neutrinos produced in a long duration 
flare from microquasars similar to, e.g., XTE J1118+480, Cyg X-3 or GRO 
J 1655-40, should produce a few neutrino events in a 1 km 2 detector. In 
the following, we report the neutrino fluxes given by Distefano et al. (2002, 
private communication): 

Gear)" < E - < 5 10S GeV » 
(lib) ' (E " > 5 103 GeV) (20) 

being K = 10 4 /1.6 a constant needed to adapt Distefano et al. flux units (in 
erg cm~ 2 s _1 ) to (GeV -1 m~ 2 s -1 ). The normalization factor (4 + 51n20) is 
due to the presence of a knee in the spectrum. F v is a number characterizing 
the specified source, and is calculated in the model. Some selected sources 
are listed in Tab.|3| along with the expected v fluxes. 

In the case of microquasars containing close massive stars of OB type, 
the relativistic hadrons accelerated in the jet can also interact with the dense 
wind of the massive star. In their model, Romero et al. (2003) consider as 
likely parameters of such a binary system the following ones: mass loss rate 



dN v _ K 
~d~E~ v ~ (4 + 5 In 20) 

5000 K 
~ (4 + 5 In 20) 
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Table 3 List of microquasars considered in Distefano et al. (2002), including 
expected neutrino fluxes in units of erg s" 1 cm" 2 . 



Source 


F 


1/ 








XTEJ1748-288 


3 


07 


X 


10" 


10 


CygnusX-3 


4 


02 


X 


10" 


9 


LS5039 


1 


69 


X 


10" 


12 


GROJ1655-40 


7 


37 


X 


10" 


10 


GRS1915+105 


2 


10 


X 


10" 


10 


CircinusX-1 


1 


22 


X 


10" 


10 


XTEJ1550-564 


2 


00 


X 


10" 


11 


V4641Sgr(l) 


2 


25 


X 


10" 


10 


V4641Sgr(2) 


3 


25 


X 


10" 


8 


ScorpiusX-1 


6 


48 


X 


10" 


12 


SS433 


1 


72 


X 


10" 


9 


GS1354-64 


1 


88 


X 


10" 


11 


GX339-4 


1 


26 


X 


10" 


9 


CygnusX-1 


1 


88 


X 


10" 


11 


GROJ0422+32 


2 


51 


X 


10" 


10 



of massive star ~ 10~ 5 M Q yr -1 , wind velocity 2500 km s -1 , radius of 
the massive star 35 Rq, mass of the black hole 10 Mq, accretion rate onto 
the black hole 10 -8 M Q yr _1 . They conclude that, if the spectrum of the 
accelerated protons is of power law type, with spectral index close to 2, and 
it extends up to 100 TeV, the fluxes of neutrinos produced in collisions of 
jet accelerated hadrons with the matter of the wind might be a factor of 3 
above the atmospheric neutrino background in the energy range 1-10 TeV, 
if the efficiency of energy conversion of the accretion disk to the relativistic 
protons is equal to 10~ 2 . 

Simultaneous observations of the X-ray and neutrino emission from spe- 
cific microquasars can help to distinguish between these two models or at 
least estimate their relative importance. 



9 Magnetars 

Magnetars are isolated neutron stars with surface dipolc magnetic fields 
much higher than in ordinary pulsar, typically ~ 10 15 G. Such objects mani- 
fest themselves in the form of soft gamma-ray repeaters (SGRs) and anoma- 
lous X-ray pulsar (AXPs). The prime distinction between a pulsar and such 
a magnetar is that the X-ray and particle emissions from the magnetar are 
powered not by rotation, as in common pulsars, but also by the decaying 
magnetic field. This involves both internal heating and seismic activity that 
shakes the magnetosphere and accelerate particles. This gradual release of 
energy is punctuated by intense outbursts that are most plausibly triggered 
by a sudden fracture of the neutron star's rigid crust. 
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In a recent paper, Zhang et al. (2003) propose that young fast-rotating 
magnetars are likely to be significant TeV neutrino sources, and therefore 
interesting targets for the planned km 3 neutrino detector. In their model, 
Zhang et al. (2003) argue that the maximum potential drop through the 
magnetar magnetosphere might be high enough to accelerate protons to 
energies above the photomeson threshold. These protons will then interact 
with the thermal radiation from the neutron star surface, thus producing 
pions and neutrinos. In order to accelerate protons to the required energies, 
the parameters of the magnetar have to fulfill the following condition 

P < (2.4 - 6.8 s) Bl%Rl ,2 {l - cosVr) 1 / 2 (21) 

where P, B p , and R are the rotation period, surface (dipolar) magnetic field, 
and stellar radius, and 9 P1 is the interaction angle between the photon and 
the proton. The range of periods in the right-hand side of eq. (|21(l defines a 
"neutrino death valley" in the magnetar P — P (or P — B p ) space. Photome- 
son interactions and neutrino emission cease when the magnetar crosses this 
valley from left to right during its evolution (see Fig.l of Zhang et al. pa- 
per). Four magnetars are found to be close to the valley, which means that 
under favorable conditions, they are high energy neutrino emitters. For an 
on-beam observer, and assuming that the neutrino luminosity is beamed 
into a solid angle AQ V around the polar axis, Zhang et al. find that the 
neutrino flux at Earth is 




(22) 

where D is the distance to the source. Moreover f coo i is the pion cool- 
ing factor, kT ~ (0.4-0.6) keV is the observed blackbody temperature for 
SGR/AXP quiescent emission, and r\ p parametrizes the uncertainty in the 
utilization of the polar cap unipolar potential. Zhang et al. have found that 
for on-beam detections, SGR1900+14 and 1E1048-5937 have substantial 
neutrino fluxes, making them interesting targets for the planned large area 
Cherenkov detectors. Unfortunately, event rates strongly depend on this 
beaming angle. In fact, a smaller Afl v increases event rates, as neutrinos 
flux depends on (\j Afl l/ ) 1 but decreases the probability of on-beam detec- 
tion. Moreover, we notice that the neutrino flux spectral index is the same 
for all sources, while the normalization factor changes as it depends on the 
distance, the rotation period and the surface magnetic field of the source. 

The estimated neutrino event rate from one of the more likely neutrino 
source, SGR1900+14 should produce ~ (1.5 - AQ y ) events in 1 yr 

and in a km 2 detector, if the protons are accelerated with the efficiency cor- 
responding to the Goldreich & Julian (1969) density. Such a rate is above 
the sensitivity of a 1 km 2 detector, provided that the observer is located 
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within the cone of neutrino emission. There might be some problems with 
acceleration of hadrons to the required energies in this model. The authors 
estimate the mean free paths for pion and pair production in hadron- 
photon collisions for typical parameters of thermal radiation field expected 
in magnetars and conclude that the electric field in the acceleration region is 
not saturated by created pair plasma. However, it is assumed that hadrons 
are injected with the Goldrcich and Julian (GJ, 1969) density and the num- 
ber of pairs produced by hadrons per one pion should be of the order 
of ten (Zhang et al. 2003), i.e. ten times above GJ density. These pairs are 
created in strong radiation and electric fields so they should initiate very 
efficient cascades since the mean free paths for the ICS and subsequent 
pair production are approximately one thousand times shorter than for pion 
production. Therefore the number of secondary pairs in such dense radia- 
tion and strong electric fields can be enormous (see e.g. Bednarek & Kirk 
1995). 

Another problem of that model is that, in the above discussions about 
magnetar gap acceleration, a dipole configuration is assumed, whilst a mag- 
netar magnetosphere is certainly not dipolar. On this point, Thompson et 
al. (2002) argue that the SGR/AXP phenomenology can be consistent with 
the hypothesis that the magnetar magnetosphere is instead globally twisted. 

10 Summary and conclusions 

In general we can distinguish transient, persistent, and periodic neutrino 
sources. Some of these sources should produce strongly collimated neutrino 
beams, e.g. supernova models (Protheroe et al. 1998), binary system models 
(Anchordoqui et al. 2003b, Bartosik et al. 2003), or magnetar model (Zhang 
et al. 2003). Therefore, a specific geometry is required in order to have a 
chance to detect them. 

The neutrino fluxes expected from transient sources such as the super- 
novae can produce from a hundred up to a few thousands of neutrino events 
in a 1 km 2 detector during the first year after supernova explosion within our 
Galaxy, provided that the accelerated hadrons take a significant part of the 
available supernova energy (Waxman & Locb 2001, Beall & Bednarek 2002, 
Nagataki 2004). However, a more realistic estimation of the hadron accel- 
eration efficiency, which is based on the condition that the observed cosmic 
ray flux in the Galaxy should not be exceeded (Blasi et al. 2000, Beall & 
Bednarek 2002), produce significantly reduced neutrino event rates in a 1 
km 2 detector. A few events can be expected in the case of pulsar formation 
with typical values of the surface magnetic fields of radio pulsars, but with 
initial periods of the order of a few milliseconds. If a pulsar with extremely 
high surface magnetic field is born in a supernova explosion, a so called 
magnetar, then the estimated event rates might be higher. However, in this 
case the pulsar can additionally accelerate the supernova envelope, dropping 
faster the column density of matter for relativistic hadrons and consequently 
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lowering the expected neutrino event rates. This effect has not been taken 
into account in the calculations presented by Beall & Bednarek (2002) and 
Nagataki (2004). In general, detection of neutrinos from the early stage of 
supernovae in a reasonable time in not very promising due to the low rate 
of supernova explosion in the Galaxy ~ 10~ 2 yr -1 . Estimated event rates 
from Galactic supernovae do not give also much hope of detecting such 
supernovae in the nearby galaxies. For example, the expected event rates 
from the supernova SN1987A in LMC, at distance of ~ 5 times larger than 
a typical Galactic supernova, are below or at the level of sensitivity of a 1 
km 2 detector in most of the considered realistic models. 

The microquasar and magnctar models (Levinson & Waxman 2001, 
Zhang ct al. 2003) predict efficient production of neutrinos during the out- 
bursts observed in the radio and/or X-ray bands. The predicted strong cor- 
relation between the expected neutrino fluxes and this low energy emission 
should eventually help in extracting the neutrino signal from the atmo- 
spheric neutrino background. This will require long multi-wavelength cam- 
paigns since the expected detection of the neutrino signal from , e.g. GRS 
1915+105, will need the observation of about 30 outbursts. 

More promising detection of neutrinos is expected from persistent sources 
such as the PWNe and the shell-type SNRs, provided that a significant part 
of their TeV 7-ray emission have hadronic origin. The expected neutrino 
event rates from these sources are relatively low in a 1 km 2 detector, of the 
order of a few events yr _1 (Alvarez-Muniz & Halzen 2002, Bednarek 2003a, 
Amato et al. 2003), but we know exactly where to look for them. Therefore, 
even a relatively weak signal from the direction of the Crab Nebula or TeV 
7-ray shell-type SNRs will give important constraints on hadronic processes. 
Another likely persistent neutrino source is expected in the Galactic Cen- 
tre region. If excess of CRs from the Galactic Centre at ~ 10 18 eV, found 
in the AGASA and SUGAR data, is confirmed (Hayashida et al. 1999, 
Bcllido ct al. 2001), then it should be accompanied by a neutrino source 
able to produce up to several neutrino events in a 1 km 2 detector per year 
(Bednarek 2002, Anchordoqui et al. 2003a, Biermann et al. 2004). In this 
case, neutrino fluxes would be correlated to the measured CR flux. 

Also the models involving hadron acceleration inside massive binary sys- 
tems are expected to be persistent and periodic (with the orbital period) 
neutrino sources. Bartosik et al. (2003) consider in more detail the old sce- 
nario in which energetic pulsar is responsible for acceleration of hadrons 
which interact with the massive companion matter. The model predicts 
that detectable neutrino fluxes can be observed, also from non-eclipsing bi- 
nary systems, due to the propagation of relativistic hadrons in the magnetic 
field of the massive star, which allows them to fall onto the star at large 
angles to the plane of the binary system. Anchordoqui et al. (2003b) argue 
that the relativistic hadrons impinging on the surface of the accretion disk 
around the neutron star, which is companion of the massive star, can pro- 
duce a neutrino signal above the atmospheric neutrino background in a 1 
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km 2 detector near the periastron for very dense accretion disks, e.g. in the 
case of A0535+26. 

Models for astrophysical neutrino sources in our Galaxy usually, very 
optimistically, assume efficient transfer of the energy generated inside the 
source into relativistic hadrons. It looks that the acceleration efficiency is 
going to be constrained for some sources by the upper limits derived on the 
basis of AMANDA observations (Ahrcns et al. 2003, Hauschildt et al. 2003). 
For example, the microquasar model (Distefano et al. 2002) prediction is 
about a factor of 2 below the AMANDA-II 90% c.l. upper limit for the 
persistent source SS433 (Ahrens et al. 2004a). 

In general, models which postulate the production of neutrinos in hadron 
collisions with the radiation field already inside the acceleration region (e.g. 
the early stage of the model by Beall & Bednarek (2002), the microquasar 
model by Levinson & Waxman (2001), or the magnctar model by Zhang et 
al. (2003)) may have problems with efficient conversion of energy from the 
acceleration mechanism into relativistic hadrons. This is due to the fact that 
the threshold for e ± pair production is about two orders of magnitude below 
the threshold for pion production and the cross section for e ± production 
is significantly larger than that for pion production in the case of hadron- 
photon collisions. Therefore, pair production is more efficient than pion 
production by protons, with the possibility to have pair domination 
of accelerated plasma. Such e ± pair domination of the plasma during the 
acceleration process in the electric field (induced in reconnection regions) 
has been discussed by Bednarek & Kirk (1995) and in the case of shock 
accelerated particles by Mastichiadis & Kirk (1995). As a result, most of 
the energy from the acceleration mechanism is transferred to leptons and 
it is radiated in the form of electromagnetic radiation but not of neutrinos. 
In the case of models which postulate production of neutrinos in hadron 
collisions with matter the situation is more promising. The cross section for 
e ± pair production in collisions of two hadrons is significantly lower (even 
for nuclei as heavy as iron) than the cross section for pion production. 
Therefore, the domination in the acceleration region by the e ± plasma is 
not possible, provided that other processes, e.g. leptonic, do not contribute 
significantly. Moreover, in these models the region of hadron acceleration is 
usually different than the region of their interaction with matter. 
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Table 4 Brief summary of the models described in this review paper. The source type and distance are indicated, the process that 
produces neutrinos, the duration of the neutrino emission (D) or/and the source age (A), the relevant energy region for producing the 
given event rates in a km 2 detector per year and the neutrino flux spectrum. 
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